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Abstract. The analysis of the stability properties of as- 
trophysical jets against Kelvin-Helmholtz (or shear-layer) 
instabilities plays a basic role in the understanding the ori- 
gin and physical characteristics of these objects. Numeri- 
cal simulations by Bodo et al. (1998) have shown that the 
three-dimensional non-linear evolution of KH instabilities 
in supersonic jets is substantially faster than in the two- 
dimensional case, leading to a cascade of modes towards 
smaller scales and a very effective mixing and momen- 
tum transfer to the ambient medium. On the other hand, 
Rossi et al. (1997) and Micono et al. (1998) found, in two 
dimensions, that radiative losses tend to reduce and de- 
lay mixing effects and momentum transfer to the ambient 
medium. In this paper, as a logical next step, we investi- 
gate the effects of radiative losses on the stability of 3D 
supersonic jets, assuming that the internal jet density is 
initially lower, equal and higher than the ambient medium, 
respectively. We find that light and equal-density radiative 
jets evolve in a qualitatively similar fashion with respect 
to the corresponding adiabatic ones. Conversely, we note 
substantial differences in the evolution of heavy jets: they 
remain more collimated and do not spread out, while the 
momentum gained by the ambient medium stays within 
~ 5 jet radii. 
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1. Introduction 

Jets that originate from young T Tauri stars in giant 
molecular clouds and from Active Galactic Nuclei propa- 
gate, in a collimated fashion, up to distances ~ 100— 1000 
times their average radii. The question of their endurance 
against Kelvin-Helmholtz instabilities is a crucial one for 
understanding the phenomenology of these objects. This 
issue has been addressed by many authors first in the lin- 
ear limit (see Birkinshaw 1991 for a review) and, in recent 
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years, through the use of 2D and 3D numerical simulations 
(see Ferrari 1998 for a review). In particular, recently, 
Bodo et al. (1998) have examined the differences between 
two-dimensional and three-dimensional jet evolution, 
showing that three-dimensional jets undergo faster evolu- 
tion, enhanced mixing and larger jet momentum depletion. 
The temporal duration of the shock-dominated phase of 
the Kelvin-Helmholtz instability evolution is strongly re- 
duced, shocks that form are weaker, and the direct asso- 
ciation between structures observed in astrophysical jets 
and shocks originating as a consequence of the nonlinear 
development of the Kelvin-Helmholtz instability is more 
problematic. Moreover, the rapid disruption that three- 
dimensional jets undergo raises the question as to how 
jets can survive for the long scale lengths shown by ob- 
servations. Although relativistic bulk velocities may be a 
stabilizing factor in extragalactic jets, the stability prob- 
lem remains crucial for the nonrelativistic jets associated 
with Young Stellar Objects (YSOs), and radiation might 
be an important ingredient in this respect. 

The interest in radiation effects on the development of 
Kelvin-Helmholtz instabilities in supersonic jets was mo- 
tivated by the observation of line-emitting knots in YSO 
jets. Linear analyses of the stability of radiative jets by 
Massaglia et al. (1992) and Hardee & Stone (1997), showed 
that radiation losses typically reduce the growth rates of 
perturbations. Extending the analysis to nonlinear ampli- 
tudes, Rossi et al. (1997) and Micono et al. (1998) found 
that radiative cooling has noticeable effects on the spatial 
and temporal evolution of two-dimensional jets; in par- 
ticular, the action of cooling leads to a longer duration of 
the initial linear stage and of the shock-dominated stage 
of the instability; shocks that form appear weaker and the 
post-shock temperatures smaller; the importance of the 
mixing phase is reduced and in some case, depending on 
the jet density, it is totally absent. Studying this prob- 
lem, Downes & Ray (1998) found a somewhat contrary 
result with respect to Rossi et al. (1997): from their tem- 
poral simulations they concluded that cooling enhances 
the momentum transfer to the ambient medium, increas- 
ing the level of mixing. They suggested that the differences 
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in the two studies might be due to the different geometries 
adopted: Cartesian in their case, and cylindrical in Rossi 
et al. (1997). 

The goal of the present work is to study the effect of 
radiative losses on the evolution of the Kelvin-Helmholtz 
instability in a three-dimensional jet, with particular fo- 
cus on the issues of mixing and momentum transfer in a 
system which is not subject to any superimposed symme- 
try. Moreover, with explicit reference to jets from YSOs, 
we aim to verify whether the action of cooling is strong 
enough to prevent a fast disruption of three-dimensional 
jets, which must survive for the length and time scales that 
are observed. Specifically, we investigate the behavior of 
dense, light and equidense jets, and analyze the effects of 
increasing the amount of energy lost through radiation by 
changing the value of a radiative control parameter which 
is directly related to the ratio between the cooling and 
dynamical time scales. 

The plan of the paper is the following: in Section 2 
we examine the physical problem and the equations used; 
the numerical scheme is described in Section 3; Section 
4 summarizes the results of our previous studies on the 
stability of adiabatic jets; the results of our simulations 
are discussed in Section 5. Finally, a summary and the 
conclusions are given in Section 6. 

2. The physical problem 

We study the evolution of a 3D fluid jet, subject to radia- 
tive losses; the equations regulating our system are thus 
the ideal fluid equations for mass, momentum and energy 
conservation: 
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where p, p, and v are pressure, density and velocity, 7 is 
the ratio of the specific heats, and C and Q are respectively 
the energy loss term and the heating term (assumed here 
constant, and equal to C taken at t — 0); the right hand 
side of Eq. || vanishes in a purely adiabatic case. 

We also follow the evolution of the number fraction of 
the neutral hydrogen atoms / n : 



df 

+ (V • V)/ n - n e [-( Ci + Cr )/n + Pr] , 



(4) 



where n c is the electron density, and Cj and c r are the 
ionization and recombination coefficients, which depend 
on the temperature T and on the electron density n e in 
the following way: 



(Rossi et al. 1997). 

Finally, we trace the jet material by means of a scalar 
field T whose evolution is followed by integrating the 
scalar advection equation: 



dT 

- + (v.V)T = 



(7) 



Cooling is treated here as in our two previous papers on 
2D cooling jets, and the reader is referred to Rossi et al. 
(1997) for a detailed description. The assumed losses are 
a reasonable approximation for temperatures T< 35,000 
K and shock velocities up to about 80 km s _1 ; both con- 
ditions are generally fulfilled in YSO jets. We recall, in 
particular, that we assume here that our jet is composed 
of atomic gas with solar abundances, and we neglect the 
formation of molecules and molecular emission. This as- 
sumption is consistent with our choice for the initial jet 
temperature, which is well above the dissociation temper- 
ature for the hydrogen molecule. 

3. The numerical setup 

3.1. The integration method 

The hydrodynamical equations are integrated numeri- 
cally using a three-dimensional version of the Piecewise 
Parabolic Method code (Colella & Woodward 1984). The 
code was parallelized through the MPI package. Multidi- 
mensionality, as well as the inclusion of radiative losses, 
are achieved through operator splitting: the PPM code 
solves the Eqs. [l] - ^) and |?]) in a conservative form; where 
non-conservative terms are present, the physical quantities 
are updated at the end of the main time step. In addition, 
the integration time step is constrained by the Courant 
condition, and is further corrected so that the tempera- 
ture does not vary by more than 10% in a single step. 

3.2. Integration domain, initial and boundary conditions 

Our 3D domain {(0,D) x {-R,R) x (-R,R)}, with D = 
lOna and R — 6.7a where a is the jet radius, is covered 
by a 256 x 256 x 256 grid, and is described by a Cartesian 
coordinate system (x, y, z). 

The initial flow structure is a cylindrical jet, with its 
symmetry axis lying along the x-direction and defined 
by (y = 0, 2 = 0). The initial jet velocity is along the 
x-direction. The interface between the jet and the sur- 
rounding material is not a vortex-sheet-like interface, but 
a smoothly varying velocity shear layer. The form of the 
initial jet velocity profile is thus: 



c; = 1.08- 10~ 8 • 13.6" 2 
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where Vq is the velocity on the jet axis, and m is a param- 
eter controlling the interface layer width; in these calcula- 
tions we set m = 8 which implies a shear layer thickness 
of approximately 0.4 jet radii. 

The initial structure of the density in our domain is 
defined by: 



po{y,z) 

POJet 



(v — 1) sech 



vV 



(9) 



direction is instead uniform, and the boundary conditions 
are periodic. This choice implies a temporal approach to 
the study of the instability evolution, as opposed to a spa- 
tial approach where inflow and outflow conditions are ap- 
plied at the left and right boundary respectively, along 
the jet flow direction. We refer to Rossi et al. (1997) and 
Bodo et al. (1998) for a more detailed discussion of bound- 
ary conditions, grid settings and comparison between the 
temporal and spatial approaches. 



where v is the ratio of the density at r = \/ y 2 + z 2 = oo, 3.3. Scaling and parameters 



to that on the jet axis (y — z — 0) at the initial time. 

We assume that the jet is initially in ionization equi- 
librium, and in pressure equilibrium with its surroundings 
(with an initially uniform imposed pressure distribution). 

Finally, we impose the following initial conditions for 
the passive tracer T: 



T 



r < 1 
r > 1. 



(10) 



so that we are able to follow the jet material during its 
subsequent evolution. 

This initial configuration is then perturbed at t = 0; 
we perturb the transverse velocities V y and V z , so that a 
wide range of modes can be excited. In order to mainly 
excite the helical mode, we choose the following functional 
form for the transverse velocities: 



V y (x, y, z) = V y fi sech 



Vv 2 + 



i " 
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V z (x,y,z) = V zfi sech 



i " 

x - ^cos(jfc x + <f>j), 
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3 = 1 



where V V: o — V z _o = 0.05Vb is the amplitude of the ini- 
tial perturbation, {</>,,} are the phase shifts of the various 
Fourier components, n — 8 and ko = 2tt/D. 

As regards the behavior at the boundaries, we adopt 
free outflow boundary conditions at the upper and lower 
boundaries in the y and z directions, by setting the gradi- 
ent of every variable to zero. The central region of the do- 
main, —3a < y < 3a, —3a < z < 3a, is covered uniformly 
by 150 grid points in the radial direction (with the jet it- 
self initially occupying 50 grid points at t = 0); at larger 
distances from the jet axis, both the y and the z mesh sizes 
increase accordingly to the scaling laws Ay J+1 = 1.02 A?/j 
and Azfc +1 = 1.02A^ respectively. The grid along the x 



Our calculations are carried out in non-dimensional form. 
Distances are thus measured in units of the jet radius a, 
velocities are scaled to the isothermal sound speed c s , 
evaluated on the jet axis and at the initial time t = 0. 
Accordingly, time is then measured in units of the sound 
crossing time t cr . 

The dynamical control parameters for our problem are 
the Mach number M(= Vo/(c sy /j)) and the density ratio 
v = Poc/Pjct- The treatment of cooling introduces two 
more parameters in our calculations, namely the initial 
temperature on the jet axis To and the ratio r ra d between 
radiative cooling time t ra d = _p/[(7 — 1)£] and the sound 
crossing time. Fixing r ra( j is analogous to fixing the prod- 
uct of the jet particle density and the jet radius n^a, which 
was selected in our calculations since it has more direct 
observational meaning (Rossi et al. 1997). 

In our simulations, we fixed the Mach number M = 10 
and we selected three values for the ratio v between the 
external and the jet material, i.e. v = 0.1 (overdense jet), 
v = 1 (isodense jet) and v — 10 (underdense jet). The 
initial temperature on the jet axis was always To = 10, 000 
K, and for the product of the jet particle density and the 
jet radius we selected the value n^a — 5 x 10 17 cm -2 , which 
is consistent, for example, with a jet radius a — 5x 10 15 cm 
and a jet particle density no = 100cm -3 ; however, we 
recall that our results are valid for all pairs of values of no 
and a for which their product is constant. 

For the equidense jet {y — 1) case we performed a 
simulation with the value n^a — 5 x 10 18 cm -2 in order to 
study the effect of varying the radiative cooling time on 
the instability evolution. 

In Table § we summarize the simulations discussed in 
the present paper, with the values of the relevant param- 
eters. 

Table 1. 



Simulation 


M 


V 


To 


no a 


A 


10 


0.1 


10 4 K 


5. x 10 1 ' cm" 2 


B 


10 


10. 


10 4 K 


5. x 10 17 cm -2 


C 


10 


1. 


10 4 K 


5. x 10 17 cm -2 


D 


10 


1. 


10 4 K 


5. x 10 18 cm -2 
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4. Results of previous studies 

4--1- Stages of the instability evolution 

Previous studies of the temporal (Bodo et al. 1994, 
1995, Rossi et al. 1997) and spatial (Micono et al. 1998) 
evolution of the Kelvin-Helmholtz instability in two- 
dimensional jets allowed us to distinguish four stages of 
evolution; they were identified in adiabatic and cooling 
jets in the spatial and temporal approach, in slab and 
cylindrical geometry; their features, as for example the 
onset time, the duration, the morphological and physical 
evolution of the jet during each stage, were dependent on 
the parameters adopted, although the general trend was 
common to all studied cases. We recall in the following 
the main features of each evolutionary stage: 

— Stage 1: the unstable modes excited by the perturba- 
tion grow in accord with the linear theory. In the latter 
portion of this stage, the growth of the modes leads to 
the formation of internal shocks. 

— Stage 2 : the growth of internal shocks is accompanied 
by a global deformation of the jet, which thus drives 
shocks in the external medium. These shocks can carry 
momentum and energy away from the jet, and transfer 
them to the external medium. 

— Stage 3 : as a consequence of the shock evolution, mix- 
ing between jet and external material begins to occur. 
The longitudinal momentum, which is initially con- 
centrated inside the jet radius, is spread over a much 
larger region by the spreading of the jet material. 

— Stage 4- in this final state, the jet attains a statistically 
quasi-stationary state. 

The effects of radiative losses on the instability have 
been traced by Rossi et al. (1997) and Micono et al. (1998). 
The main differences of this analysis with respect to the 
adiabatic case are: i) a longer duration of stages 1 and 2 
of the evolution; ii) a decrease of the shock strength and 
of the post-shock temperatures, and iii) a general reduc- 
tion of the importance of the mixing phase (depending on 
the values of v and hqo), with a complete absence of the 
mixing phase for high values of no a. 



4.2. 2D vs 3D behaviour 

In this section we summarize the main results obtained by 
Bodo et al. (1998) in their comparative analysis of 2D and 
3D adiabatic jets. Bodo et al. (1998) adopted a Cartesian 
geometry and a temporal approach (i.e., periodic bound- 
ary conditions in the direction of the jet axis), and super- 
imposed helical perturbations onto the initial transverse 
velocities (see also below). 

The expected differences are mainly due to two fac- 
tors: first, there are many more unstable modes in three 
dimensions than in two dimensions, and the growth rates 
of the 3D modes can predominate. Second, it is well-known 



that turbulence in two and three dimensions differs sub- 
stantially, especially as far as the mixing properties are 
concerned. Both effects - the presence of unstable high- 
wavenumber modes, and the cascade to high wave num- 
bers via nonlinear interactions - can lead to the formation 
of small-scale structures that are not observed in two di- 
mensions. 

The main features observed by Bodo et al. (1998) 
in the evolution of over-dense and under-dense three- 
dimensional adiabatic jets, compared to analogous two- 
dimensional slab jets, can be summarized as follows: 

1. Faster evolution. 

2. Absence of a clear separation of stages 2 and 3 of the 
instability (as defined in section 4.1). 

3. More rapid development of small-scale structures due 
to the growth of linearly unstable higher order (az- 
imuthal number < 2) fluting modes and, at advanced 
times, to a cascade to small scales through nonlinear 
turbulent processes. 

4. Enhanced material mixing between the jet and the am- 
bient medium; more effective and faster momentum 
transfer from jet to the ambient medium. 

5. Enhanced jet broadening. 

6. Similar evolution of the largest-scale structures. 

7. Similar asymptotic states (although reached at very 
different times). 



5. Results 

5.1. Heavy jets (v = O.lj 

In this section we present the results of the calculations 
for an adiabatic jet and a cooling jet with Mach number 
M = 10 and density ratio v = 0.1. We carried out the 
simulation for t = 12 units of sound crossing time, since 
at this time the adiabatic jet reaches the boundaries in 
the y and z directions. 

In Fig. [l] we show grey-scale images of two-dimensional 
cuts of the jet density, in the z = plane (i.e., a longitu- 
dinal cut in the (xy) plane in the middle of the domain) 
and in the x — 5tt plane (i.e. a transverse cut in the (yz) 
plane in the middle of the domain), for the radiative and 
adiabatic jets at different stages of the evolution. 

At time t = 4 the situations is similar for the adia- 
batic and cooling jets: the instability has already entered 
the acoustic phase, with the development of shock waves 
driven by the jet into the external medium, in correspon- 
dence to the major "kinks" caused by the growth of the 
fundamental helical mode excited by the imposed pertur- 
bation. Acoustic waves are driven in the external medium, 
as it can be clearly seen in the upper panels of Fig. [I]. 
Already at this early stage, a little amount of external 
material is being entrained by the jet, as a result of the 
growth of higher order fluting modes that produce short 
scale ripples on the surface of the jet. 




Actiababic time — B.O Radiativa time — S.O 




5 tO fS SO ES 30 -f+SHSJS £f G iO IG SO SS 30 -B-4-B0B4-6 



Adiabatic 



time = 8.0 



Radiative 



time = 8.0 



























0,2 








0.0 



6 io 16 20 25 20 4 8 

S B 



E 10 Iff HO £G SO -Mi^fi 

T E 



time = tg r Q 








ef 






4 i 






St 






A of 






-i. ! 










1 ...1.. .1... ..1... f 


-el 



5 tO ?5 Jfl W 30 -g-+20Z4B 




Fig. 1. Grey-scale images of the density distribution for heavy (v = 0.1) adiabatic (left column) and radiative (right 
column) jets; 2D cuts at a hxed z (z = 0) and at a fixed x (x = 5tt). 



At time t — 6 the radiative jet is still shock-dominated, 
although the disturbances induced on the jet surface have 
grown in amplitude and the amount of entrained ambient 
medium has increased. In contrast, the adiabatic jet ap- 
pears to have fully entered the mixing phase; no trace of 
the coherent large-scale deformations induced by the fun- 
damental helical Kelvin-Helmholtz mode can be discerned. 
As the evolution proceeds, the adiabatic jet appears to 
be completely mixed with the ambient medium (t = 8, 
Fig. [I]), while the jet subject to radiative losses maintains 
its collimation up to the latest time in our calculations, 
t = 12; at this time both jets have reached stage 4 of the 
instability evolution, i.e., a statistically "quasi-stationary" 
stage. The comparison between the adiabatic and radia- 



tive jets in Fig. |l| also conveys the impression that the 
adiabatic jet has become much wider than its radiative 
counterpart. To analyze this effect in more detail, we have 
computed the radial distribution of velocity averaged over 
the longitudinal and azimuthal directions and we have de- 
fined the jet radius as the radius at which this averaged 
velocity drops to one half of its maximum value (found at 
r = 0). Fig. |^c) shows the evolution of this radius as a 
function of time; we see from the figure that in fact the 
adiabatic jet widens much more than the radiative one. 

In Fig. we show a two-dimensional cut of the Lapla- 
cian of the density in the plane z = (the Laplacian of the 
density is a valuable tool to mark the presence of shock 
waves in the flow, and is widely used in gas dynamics). 
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Fig. 2. Grey-scale images of two dimensional cuts of the density Laplacian, in the plane z = 



Upon examining this figure, we can follow the same evo- 
lutionary path described above. The Laplacians at time 
t = 2 arc not shown since no significant feature can be de- 
tected: both jets are in fact still in the linear phase (stage 
1) of the evolution. When the perturbations steepen into 
shock waves, the jet enters the acoustic phase (t = 4). The 
duration of the acoustic phase is longer for the radiative 
jet, where strong shock waves survive up to t = 10; for 
the adiabatic jet, instead, shock waves are damped by the 
vigorous mixing that occurs already during the acoustic 
stage: after time t — 6, no shocks can be detected in the 
flow, and turbulent mixing dominates. 

In panel (a) of Fig. ^ we plot the tracer entropy W vs. 
time, for the adiabatic (dashed line) and radiative (solid 
line) dense jets. W measures the departure of the tracer's 
distribution from the initial form (T = 1 for the jet ma- 
terial and T = for the external material; see Bodo et al. 
(1995) for the mathematical definition of W). Up to time 
t = 2 the jet is still in the linear phase of the evolution, 
and the entropy is almost constant, since the tracer distri- 
bution does not differ substantially from the initial one. 
As time elapses, the jet goes through stages 2 and 3 of 
the evolution, during which the entropy grows, reaching a 
quasi-constant value at the end of the evolution. 

In Fig. ^|b) we plot the momentum of the radiative 
(solid line) and adiabatic (dashed line) jets, normalized to 
the initial values, as a function of time. In both cases, the 
jets start losing momentum early in time; the final momen- 
tum loss is dramatic, since at time t = 12 radiative and 
adiabatic jets have lost almost 80% of their initial momen- 
tum. Fig. |3p) shows that the amount of momentum lost by 



the radiative and the adiabatic jets is similar; the momen- 
tum loss is slightly greater for the radiative jet. This result 
seems to be in contradiction with the behavior described 
previously: if we assume that the quantity of momentum 
lost is a valid indicator of the instability effects, then we 
would expect that the radiative jet is equally or more un- 
stable than the adiabatic jet; on the other hand, we saw 
that radiative cooling acts to reduce the development of 
mixing and to conserve the collimation of the jet for longer 
times. In fact, as we have seen above, the jet widening in 
the adiabatic case is larger than in the radiative case; one 
could then expect that the amount of entrained material 
would be much larger in the adiabatic case, and corre- 
spondingly the residual jet momentum much lower. This 
apparent discrepancy can be explained if we analyze how 
momentum is lost by the jet in the two cases. 

The rate of momentum transfer is almost similar for 
the cooling and adiabatic jets, but the radiative jet accel- 
erates only the ambient medium located very near the jet 
itself; in contrast, the adiabatic jet expands and entrains 
material over a much larger radial extent, accelerating it 
until at the end of our calculations all the ambient medium 
in our computational domain has a non-zero velocity in 
the longitudinal direction; this trend appears clearly in 
Fig. ||a), where the amount of momentum acquired by the 
external medium at different times is shown as a function 
of the distance from the initial jet axis. 

The typical velocities attained by the accelerated am- 
bient medium in the two cases are similar (for example, 
the maximum velocity of the material with T < 0.1 at 
t=10 is ~ 6.8 for both cases), but the adiabatic jet acceler- 
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Fig. 4. Panel (a): ambient medium momentum as a function of the distance from the jet axis. Panel (b): ambient 
medium momentum versus the corresponding value of the density. In both panels we plot the trends referring to heavy 
radiative (solid line), and adiabatic (dashed line) jets at different times. 



ates a larger radial volume of ambient material. A similar 
amount of transferred momentum is thus achieved in the 
two cases only if in the radiative case denser ambient ma- 
terial is accelerated. That this is the case is demonstrated 
by Fig. ^b), which shows the distribution of momentum 
in the external medium versus the corresponding value of 
the density for the radiative and adiabatic jets: we char- 
acterized each volume of ambient medium by its values 
of momentum and density, divided the density range into 
density bins, and integrated the momentum over each den- 
sity bin, obtaining the distribution showed in Fig. ||b. This 
figure shows that the typical densities of the entrained ma- 
terial are much larger in the radiative case than in the adi- 
abatic case; furthermore, the total mass of the entrained 
material turns out to be of the same order in both cases. 
We can understand this if we note that the initial shocks 
driven by the adiabatic jets heat the external gas, which 
consequently expands and thus lowers its density. In the 



radiative case, cooling suppresses the expansion, leading in 
contrast to a compression of the shocked material; the mo- 
mentum acquired by the external medium resides largely 
in material whose density is up to ten times higher than 
the initial external density. In addition we must remember 
that this shock-induced compression can be very effective 
in the radiative case since strong shock waves survive for 
long times, from t = 4 to t = 10 (see Fig. ||). 

5.2. Light jets (u = 10) 

A quick glimpse at Fig. |5|, which displays two-dimensional 
grey-scale cuts of the density at different times in the z = 
and x = 57r planes, shows that radiative losses in under- 
dense jets have little effect on the instability evolution, 
which is slower with respect to the previously examined 
cases: the linear phase of the instability (stage 1) lasts 
until time t = 8 when the first shocks appear in both the 
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Fig. 5. Grey-scale images of the density distribution for light [y = 10) adiabatic (left column) and radiative (right 
column) jets; 2D cuts at a fixed z (z = 0) and at a fixed x (x = 5tt). 



adiabatic and in the cooling jets. The density distributions 
shown in Fig. || cannot show details inside the jet because 
the grey-scale saturates there. Therefore, in order to ob- 
tain better views of the jet structures, we show in Fig. 
H grey-scale two-dimensional cuts of the product of the 
tracer and the density at different times. From these fig- 
ures we can see that shocks form inside the jet and in the 
external medium, corresponding to the kinks induced by 
the fundamental helical mode, and persist for long times: 
at time t = 16, when mixing has taken place and the jet 
diameter has considerably widened, a few weak shocks can 
still be detected both in the adiabatic and cooling jets. 

The growth of the higher order fluting modes is de- 
layed in time, and mixing occurs later, mainly because of 
the inertia of the external medium. The onset of mixing 
occurs between times £ = 10 and t = 12, and is initially 
limited to the more external layers of the jet. At later 



times (e.g., time t — 18 in Fig. ||) the action of mixing is 
not disruptive, and the jet maintains its collimation. 

The transition between the evolutionary stages is well 
depicted by the tracer entropy, plotted in Fig. fj]a): en- 
tropy remains almost constant, increasing slowly for the 
first 8-10 sound crossing times, when the perturbations 
grow linearly and the first shocks form. A larger increase 
of the entropy occurs at later times, with a steeper growth 
for the adiabatic case (which reaches the quasi-stationary 
phase earlier, at t ~ 15 as compared to t ~ 18 for the ra- 
diative case). The jet momentum is plotted in Fig. f?]b) as 
a function of time; we can see that the momentum behav- 
ior is similar to the entropy behavior: the drop in the jet 
momentum appears to be steeper for the adiabatic case, 
but it stops earlier so that the total amount of momentum 
lost by the jet material is about the same for the two cases 
(~ 90% of the initial value). A more careful comparison 
between the behaviors of entropy and momentum shows 
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Fig. 6. Grey-scale images of the jet density (i.e., of the product between density and tracer) distribution for light 
(y = 10) adiabatic (upper line) and radiative (lower line) jets; 2D cuts at z = 0. 



that the phase of steeper decrease of momentum begins 
earlier than the phase of steeper increase of entropy. In this 
case we then have a distinct phase 2 (acoustic phase) of 
the evolution in which momentum is mainly lost through 
shocks. 

The largest differences between the adiabatic and the 
radiative cases can be seen in the last stage of the evolu- 
tion: looking at Figs. || and || at times £ = 16 and t = 18, 
we see that the adiabatic jet appears to be wider than the 
radiative one. This is confirmed by Fig. 0c), which shows 
the behavior of the jet radius, as defined in the previous 
subsection, as a function of time. The difference between 
the adiabatic and the radiative cases however is not so 
large as for the dense jets. This same effect can also be 
seen in the jet momentum distribution: in Fig. || we show 
the grey-scale distribution in the (yz) plane of the aver- 
age of the jet momentum over the longitudinal direction 
for the two cases; we see that the distribution is wider for 
the adiabatic jet. Since the total amount of momentum 
retained by the jet is about the same in the two cases, we 
also find that the maximum value of the jet momentum, 
located at a position corresponding to the initial jet axis, 
is smaller for the adiabatic jet (about half the maximum 
value of the cooling jet momentum). 

From the above description we can conclude that - in 
contrast to the heavy jet case - the presence of radiative 
losses does not strongly modify the overall evolution of 
the instability in a light jet. To investigate the reasons 
for this behaviour, we plotted in Fig. ^ the total radia- 
tive power (jet plus ambient medium) for heavy, light and 



equal-density jets. In each case, the peak in the emission 
takes place shortly after the onset of the acoustic stage, 
when shocks are well developed and not yet disrupted by 
mixing. 

The amount of energy lost through radiation is at least 
an order of magnitude smaller in the light jet compared 
to the dense jet. The reasons for this behavior are many; 
first, the growth of the helical mode in the dense jet drives 
larger amplitude kinks when compared to those observed 
in the light jet; in this way the shock fronts are wider in the 
v = 0.1 jet, and a larger amount of material is compressed, 
heated and ultimately cooled through radiation. In addi- 
tion, the wavelength of the most unstable mode may play 
an important role: a shorter wavelength, typical of dense 
jets (see Hardee & Norman 1988) leads to the formation 
of a higher number of shocks in a short time, increasing 
the efficiency of the cooling mechanism at shocks (com- 
pare also Figs. [I] and ||). We also note that the ambient 
medium gives a significant contribution to the subtraction 
of energy from the system, in particular when dense and 
equal-density jets are concerned; in fact, shocks are driven 
by the jet into the ambient medium, and they propagate 
in the external regions immediately surrounding the jets; 
the spatial extension of these shocks can be consistently 
larger that in the jet's body. 

5.3. Equal- density jets (v = 1 ) 

When studying equal-density radiative jets with the same 
value of n a as considered in the previous cases (n a = 
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Fig. 3. Panel (a) shows a plot of the tracer entropy for the 
radiative (solid line) and adiabatic (dashed line) jets as a 
function of time (for the analytical and physical definition 
of tracer entropy see Bodo et al. 1995); panel (b) shows a 
plot of the jet momentum for the radiative (solid line) and 
adiabatic (dashed line) jets as a function of time. Panel (c) 
provides a plot of the temporal evolution of the jet radius 
for the radiative (solid line) and adiabatic (dashed line) 
cases . All the plots are for the heavy jet case (v = 0.1). 
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Fig. 7. In panel (a) we show a plot of the tracer entropy 
for the radiative (solid line) and adiabatic (dashed line) 
jets as a function of time; in panel (b) we show a plot of the 
jet momentum for the radiative (solid line) and adiabatic 
(dashed line) jets as a function of time. In panel (c) the 
temporal evolution of the jet radius for the radiative (solid 
line) and adiabatic (dashed line) cases is plotted. All the 
plots are for the light jet case (y = 10). 
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5 x 10 cm , consistent, for example, with a jet radius 
a = 5 x 10 15 cm and a jet particle density no = 100 
cm -3 ), we find that radiation introduces little effect on 
the instability evolution, as was the case for light jets. 
This can be seen by comparing the first two columns of 
Fig. where we show two-dimensional cuts of density in 
the plane z = for the adiabatic and radiative jets. 

The smallness of the effects introduced by radiation 
could have been predicted from the analysis of Fig. [| 
which shows that the amount of energy lost by radiation 
is similar to the light jet case, and is smaller by more than 
an order of magnitude when compared to losses associated 
with the heavy jet. For this reason, we have examined an 
additional case with noa increased by an order of magni- 
tude, i.e., noa — 5 x 10 18 cm~ 2 , which is consistent, for 
example, with a jet radius a = 5 x 10 15 cm and a jet par- 
ticle density uq — 1000 cm~ 3 . This will increase radiative 
losses; for example, the amount of energy lost by radia- 
tion at time t — 8 corresponds to 0.08% of the kinetic 
energy of the jet at that same time, compared to a value 
of 0.02% obtained in the previous case. The results for this 




Fig. 8. Longitudinally averaged distribution of the jet mo- 
mentum at time t=18 for the adiabatic and radiative light 
jets. 
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Fig. 10. Grey-scale images of the density distribution for equidense (v = 1) adiabatic (left column), and radiative jets, 
with the product of the jet density and radius n^a = 5 x 10 17 cm~ 2 (center column) and n^a = 5 x 10 18 cm~ 2 (right 
column); 2D cuts at a fixed z (z = 0). 



high-cooling case are also shown in Fig. [ITj, in the third 
column. 

Looking at the general evolution, the upper panels of 
Fig. [lO] show that at time t = 6 shocks have formed, as a 
consequence of the non-linear growth of the helical mode; 
the fluting modes drive deformations of the jet surface and 
entrainment of the external medium by the jet. (Previous 
times are not represented since t = 6 is the first time 
when non-linear effects can be detected.) The evolution 
of fluting modes, which are responsible for early mixing, is 
similar in equal-density jets and in heavy jets, for both adi- 
abatic and radiative cases: however, mixing between the 
jet and the external medium starts later (at time t = 6) 
with respect to the v — 0.1 case. Mixing has an imme- 
diate and disruptive effect: external material is entrained 



and the jet diameter widens quickly; at time t = 12 the 
jet maintains its large-scale structure and collimation, but 
on small scales it appears completely mixed with the am- 
bient medium. Comparing the results for the three cases, 
we see that the radiative cases show stronger density en- 
hancements inside the jets, and that the widening in the 
radiative cases is reduced (this reduction is particularly 
evident in the high-cooling case). This is again confirmed 
by Fig. |IT|c), which shows the behavior of the jet radius, 
as defined above, as a function of time. 

As for the heavy jet case, the acoustic and mixing 
stages cannot be distinguished in the jet evolution, both 
for the adiabatic and radiative cases. The onset of mix- 
ing is contemporary with the formation of the first shocks 
in the flow, but the interaction of the two phenomena 
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Fig. 9. Power lost through radiation by heavy (solid line), 
light (dotted line) and equal-density jets (dashed line). 
Units are normalized to the peak emission in the v = 0.1 
jet. 

is different if cooling is significant: in the adiabatic and 
no = 100 jets, vigorous mixing takes place and shocks are 
destroyed after a few sound-crossing times, while in the 
cooling n — 1000 jet strong shocks survive up to the lat- 
est stages in the jet evolution, when the global collimation 
of the jet has been partially destroyed. In fact, in the adi- 
abatic jet a strong shock appears only at time t = 6, while 
no shocks survive after time t — 8 when a weak discontinu- 
ity in the flow can be traced. In the cooling no = 1000 jet, 
instead, shock waves are not destroyed up to time t = 14; 
this is approximately the time at which the cooling jet 
enters the quasi-stationary stage of the evolution, as it is 
indicated by the behavior of the tracer entropy shown in 
Fig. [ll]a). Entropy grows more steeply for the adiabatic 
case: the quasi-stationary phase is attained earlier, i.e., at 
time ~ t = 12. 

Fig. |Ti"|b) shows the variation in jet momentum as a 
function of time for the radiative and adiabatic jets. The 
momentum loss appears to occur slightly faster for the 
adiabatic and n = 100 jets, although the final loss at- 
tained at the end of the evolution is similar for the three 
cases, and corresponds to the ~ 10% of the initial jet mo- 
mentum. Although the integrated values are similar in all 
cases, the spatial distribution of the jet momentum reflects 
the differences in the jet widening noted above. 

6. Energy spectra 

We have seen that radiation losses tend to maintain jet 
collimation. To obtain a different perspective on the insta- 
bility evolution and on the differences between the various 
cases, we have analyzed the relative amounts of kinetic 
energy in longitudinal and transverse motions at different 
scales. 

When the instability grows, the kinetic energy of the 
longitudinal jet bulk motion goes into internal energy, ra- 
diation, and kinetic energy in the transverse motions. In 
the direction perpendicular to the jet axis we can have: a) 
large-scale motions, due to the growth of the helical un- 
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Fig. 11. In panel (a) we show a plot of the tracer en- 
tropy for the radiative, no = 100 (solid line), radiative, 
no = 1000 (dotted line) and adiabatic (dashed line) jets 
as a function of time; in panel (b) we show a plot of the 
jet momentum for the radiative, no = 100 (solid line), 
radiative, no = 1000 (dotted line) and adiabatic (dashed 
line) jets as a function of time. In panel (c) the tempo- 
ral evolution of the jet radius for the radiative, no = 100 
(solid line) , radiative, no = 1000 (dotted line) and adia- 
batic (dashed line) cases is plotted. All the plots are for 
the equal density jet case {y — 1). 

stable modes, and b) small-scale motions, induced by the 
Kelvin-Helmholtz fluting modes, and by a possible turbu- 
lent cascade from large to small scale motions. Our aim is 
to analyze in more detail the transfer of energy from lon- 
gitudinal to transverse motions; the key issue is to identify 
possible differences due to radiative losses. This analysis 
can also give us informations about the persistence of the 
jet as a coherent structure of longitudinal velocity. 

We perform this analysis by calculating the three-di- 
mensional Fourier spectra of the kinetic energy for longi- 
tudinal and transverse motions. We then obtain one-di- 
mcnsional spectra by integrating first over the transversal 
wavenumbers k y and k z , and second over the longitudinal 
and azimuthal wave numbers (after transforming to cylin- 
drical coordinates). We use these one-dimensional Fourier 
spectra A x (k x ) and A r (k r ) to compare the amount of en- 
ergy in the longitudinal and radial direction, at all times, 
and for all the studied cases. Initially, before the jet enters 
the non-linear phase, the energy is found in longitudinal 
motions at large scales; as the instability evolves, a portion 
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Table 2. Value of the maximum kinetic energy in the transverse motions, time at which this maximum is reached and 
time at which equipartition between longitudinal and transverse energy is attained, for all cases, and in four different 
scale intervals. Blank spaces are left where equipartition was not observed within the maximum evolutionary time 
reached in our calculations. 
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of the kinetic energy is transferred to longitudinal motions 
at smaller scales and to transverse motions. In Fig. [l^ we 
can follow this process for the particular case of under- 
dense jets. The figure shows the behavior of longitudinal 
and transverse energies, normalized to the initial jet ki- 
netic energy, as a function of time for different scale ranges 
(recall that lengths are given in units of the jet radius). 
Looking at the longitudinal energy at large scales (A > 2, 
panel a) we see that this component, as noted above, ini- 
tially represents the only form of energy present in the 
system and, as the instability evolves, this component de- 
creases monotonically. At smaller scales (panels b,c,d) the 
longitudinal energy first increases, reaches a maximum, 
and then decreases. The time at which the maximum is 
reached is about the same for all the scale ranges. A sim- 
ilar behavior is found for the transverse energy, but the 
maximum is reached at later times. The faster decrease 
observed in longitudinal energy when compared to that 
observed for the transverse energy leads at some time to 
an equipartition between the two and, from then on, mo- 
tions can be considered to be isotropic. This isotropy is 
first reached at small scales, and then at larger scales. We 
further note, by comparing the curves for the adiabatic 
case and those for the radiative case, that in this last case 
the transfer of energy to transverse motions is less efficient 
and more energy is always kept in longitudinal motions. 

The behavior described for the underdense case is 
found also for the other cases; we have then summarized 
all our results in Table pi where we report, for all our cases 
and for each scale interval, the time at which the maxi- 
mum of transverse energy is reached, the value of this 
maximum, and the time at which equipartition between 
longitudinal and transverse energy is attained. 

By comparing the table with the behavior of entropy 
and momentum shown above, we can see that the time at 
which the maximum transverse energy is reached falls in 
the steepest portion of the momentum drop and of the en- 
ergy growth, confirming the connection between the mix- 
ing process and the transfer of energy to the transverse 
motions. In the case of radiative jets, this process of en- 
ergy transfer is less efficient as demonstrated by the lower 
values of E max for these cases and, as discussed above, 
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Fig. 12. Kinetic energy for longitudinal and transverse 
motions, as a function of time for light jets. Solid and 
dotted lines represent the energy of longitudinal motions 
respectively for adiabatic and radiative jets; dot-dashed 
and dashed lines represent the energy of transversal mo- 
tions again for adiabatic and radiative jets. 



mixing and jet widening is lower for these types of jets. 
Comparing the different values of the density ratio, we 
see that E max increases slightly as v is increased, i.e., go- 
ing from overdense to underdense jets. This increase how- 
ever cannot compensate for the increase in density with 
v\ therefore we expect lower transverse velocities for the 
high v cases, and therefore a smaller jet widening, as it is 
indeed observed in our simulations. 

The table also shows that near the end of our simula- 
tions for the v = 1 and v — 10 cases, we reach equiparti- 
tion between the longitudinal and transverse energies also 
at scales larger than the jet radius, while in the overdense 
case equipartition is found only at small scales. We can 
therefore ask whether, in these cases, the jet can still be 
identified as a coherent velocity structure. To answer this 
question, we have computed the autocorrelation function 
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for the longitudinal velocity, defined as 

F{h) = I v x (x,y, z)v x (x + h,y, z)dxdydz 



where the integral is performed over the whole domain of 
integration. We have then defined a correlation length as 



L 



F(h)/F(0)dh 



and have plotted in Fig. [13| l c as a function of time for 
all the different cases. The three panels in the figure re- 
fer respectively to the overdense case (upper panel), the 
equal density case (middle panel) and the underdense case 
(lower panel); the solid curves are for adiabatic cases and 
the dashed curves are for radiative cases. 



From Figure [13| we see that, as was also suggested by 
the Fourier analysis described above, the overdense case 
maintained an higher coherence than the other two cases, 
especially compared to the low density case in which l c 
drops to ~ 0.4 at the end of the simulation. Radiation acts 
in different directions for the different values of v: in the 
overdense case radiation tends to give a more inhomoge- 
neous structure both in density (as we have already seen) 
and in velocity and this is reflected by the lower value of 
Z c , in the equal density case the effects of radiation are 
very small, while in the underdense case radiation tends 
to keep a more coherent structure with an higher value of 

lc 

7. Summary and Conclusions 

We studied the evolution of the Kelvin-Helmholtz insta- 
bility in three-dimensional jets, solving the hydrodynami- 
cal equations through a finite-differences numerical code, 
adopting a temporal approach and varying the external- 
to-ambient medium density; we thus considered dense, un- 
derdense and iso-dense jets. The main results of our cal- 
culations can be summarized as follows: 

Heavy jet 

Dense jets undergo fast evolution. Both helical and 
fluting modes grow non-linearly at the earliest stages 
of the evolution, giving rise to large scale deformations 
in the structure of the jet and to entrainment of ambi- 
ent material. If cooling is present, the duration of the 
acoustic phase of the evolution is increased, and shocks 
survive in the flow for longer times. Furthermore, the 
strength of the fluting modes is reduced; as a result, 
mixing with the ambient medium is less effective and 
the jet retains its collimation for longer times. The mo- 
mentum loss in the adiabatic and radiative cases is sim- 
ilar, but it occurs in different ways: in both cases, the 
momentum is transferred away from the jet axis. How- 
ever, the adiabatic jet loses momentum by expansion 
and entrainment of ambient medium far from the jet 
axis, material which is heated and expands by shocks 
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Fig. 13. Longitudinal velocity autocorrelation length as a 
function of time, the autocorrelation length is in units of 
the size of the domain. The upper panel refers to v = 0.1, 
the middle panel refers to v = 1 and the lower panel refers 
to v = 10. The solid curves are for the adiabatic cases, 
while the dashed curves are for the radiative cases. 



that are driven into the ambient matter from the jet; 
in contrast, the radiative jet transfers its momentum 
to ambient material which is located near the jet sur- 
face, material which has been shocked and cooled, and 
thus condensed relative to its initial state. 
Light jet 

For underdense jets the evolution develops over 20 
sound crossing times, and thus is slower than the pre- 
vious case. Mixing takes place later with respect to the 
onset of the acoustic stage, and its action is damped 
by the inertia of the external medium. The final jet 
widening is caused by the growth in amplitude of the 
fundamental helical mode, which is the main source of 
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shock formation and momentum loss. The amount of 
energy lost through radiation at a fixed time is smaller 
than in the previous case, but cooling has nonetheless 
a small effect on the instability evolution, again reduc- 
ing the amplitude of the jet widening. 
Equidense jet 

Iso-dense jets evolve similarly to dense jets, although 
the onset of the various stages is delayed in time. The 
amount of energy lost through radiation is smaller, 
and the short duration of the shock-dominated stage 
suppresses the effects of cooling on the overall insta- 
bility evolution. To demonstrate the effects of radia- 
tive losses, we ran a simulation in which we increased 
the jet particle density (which is one of the radiative 
control parameters). The main effects of cooling arc a 
longer duration of shocks in the flow, damping of the 
helical unstable mode (and thus of the amplitude of the 
jet deformations in the transverse direction) and of the 
fluting modes, which implies that mixing is limited to 
the external layers of the jet. 

These results indicate that radiative losses indeed 
cause jet stabilization, an effect which is more evident in 
the overdense case and less evident in the other cases; they 
are thus consistent with the results of previous 2D stud- 
ies by Rossi et al. (1997) and Micono ct al. (1998), even 
though the effect is not large as in the cases studied by 
these authors. This difference can be explained as a result 
of the assumptions of (i) cylindrical geometry with axial 
symmetry and (ii) axi-symmetric initial perturbations for 
the 2D cases. The evolution of the instability in two di- 
mensions is then dominated by strong shocks on the jet 
axis: the amount of energy lost through radiation in these 
strong compressions is high, and the stabilizing effects of 
cooling are enhanced. In our 3D study, the strength of 
the shock waves that form is smaller, both because of ge- 
ometry and because of small-scale motions driven by the 
growth of the Kelvin- Hclmholtz fluting modes (the latter 
are neglected in two dimensions). In this way, the effects of 
cooling are weaker, and their importance depends on the 
density ratio between the jet and the ambient medium. 
However, our results can confirm that cooling is efficient in 
reducing the disruptive action of the instability, in main- 
taining the jet collimation for longer times, and in extend- 
ing the temporal duration of the stage where shocks are 
present in the flow. 
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